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: Abstract — Azulene-doped naphthalene was directionally solidified during the vertical Bridgman- 
j Stockbarger technique. Doping homogeneity and convection were determined as a function of the 
temperature profile in the furnace and the freezing rate. Convection velocities were two orders of 
magnitude lower when the temperature increased with height. Rarely was the convection pattern 
axisymmetric, even though the temperature varied less than 0.1 K. around the circumference of the growth 
ampoule. Correspondingly the cross sectional variation in azulene concentration tended to by asymmetric, 
especially when the temperature increased with height. This cross sectional variation changed dramatically 
I along the ingot, reflecting changes in convection presumably due to the decreasing height of the melt. 
" Although there was large scatter and irreproducibility in the cross sectional variation in doping, this 
variation tended to be least when the growth rate was low and the convection was vigorous. It is expected 
that compositional variations would also be small at high growth rates with weak convection and flat 
- interfaces, although this was not investigated in the present experiments. Neither rotation of the ampoule 
I nor deliberate introduction of thermal asymmetries during solidification had a significant influence on 
cross sectional variations in doping. It is predicted that slow directional solidification under microgravity 
| conditions could produce greater inhomogeneilies than on Earth. Combined use of microgravity and 
1 magnetic fields would be required to achieve homogeneity when it is necessary to freeze slowly in order 
to avoid constitutional supercooling. 



1. INTRODUCTION 

The Bridgman technique is widely used for direc- 
tional solidification of metals, semiconductors, inor- 
ganic salts and organic compounds. In the 
Bridgman-Stockbarger manifestation, an ampoule 
containing the material to be solidified is slowly 
translated from a furnace held at a temperature above 
the melting point of the material into a chamber held 
below the melting point. In the gradient freeze mode, 
the ampoule is held fixed in a furnace which has a 
temperature decreasing continuously from one end to 
the other. Solidification is caused by slowly lowering 
the temperature of the entire furnace. In the heat 
exchange method, coolant is applied to one end of the 
ampoule. Sometimes at the same time the furnace 
temperature is slowly decreased. 

When Bridgman crystal growth is done with the 
ampoule held vertically, it is normally desired that the 
solidification proceed upwards. Thus the temperature 
at the top of the apparatus is higher than at the 
bottom of the apparatus. For many years the crystal 
growth community referred to this arrangement as 


tPaper IAF-90-358 presented at the 41st Congress of the 
International Astronautical Federation , Dresden, 
Germany, 6-12 October 1990. 

^Present address: II-VI Inc., Saxonburg, PA 16506. U.S.A. 

357 


“thermally stable” and many persons believed buoy- 
ancy driven convection was absent. However interest 
in crystal growth in space led to a closer examination 
of this assumption. Carlson [1] at Clarkson Univer- 
sity and Brown [2] at MIT performed theoretical 
studies of convection in the Bridgman-Stockbarger 
technique with a step function change in temperature 
from the upper chamber to the lower chamber. They 
found that radial temperature gradients were in- 
evitably present in the melt and could lead to vigor- 
ous convection. Initially these results were greeted 
with skepticism by many in the crystal growth com- 
munity. Now it is widely accepted that significant 
convection is always present unless the operation is 
carried out in space or with a strong magnetic field 
applied to an electrically conducting melt. 

By now many theoretical studies have been per- 
formed on buoyancy driven convection in the vertical 
Bridgman-Stockbarger technique, examining the 
influence of gravitation level, interface curvature, 
insulation thickness between the two furnaces, am- 
poule properties, density variations caused by vari- 
ations in melt composition, varying temperatures in 
the heating and cooling chambers, etc. (e.g. [1-7]). 
Until recently [7] it was assumed that the convection 
was axisymmetric and steady. With one exception [5], 
theoretical predictions for convection have not been 
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compared with experiment because the theoretical 
work was all done for metals and semiconductors 
opaque to light. (In Ref. [5], convection was observed 
in an organic compound. Measured temperatures 
and interface shapes were used to calculate the 
convection.) 

Potts [8,9] performed experiments on an organic 
compound (high Prandtl number) heated and cooled 
by constant temperature baths separated by a dia- 
phragm. The convection in the melt was neither 
axisymmetric nor steady, although it was very weak. 
While temperature variations could not be detected in 
the top heating bath, it was significant that the 
convection always rose on the side of the ampoule 
nearest the heater submerged in the circulating bath. 

In order to see if either non-uniform heating or 
misalignment from vertical is necessary to cause the 
convection to be asymmetric, Neugebauer con- 
structed an apparatus with vertical heating wires 
embedded in plastic [10,1 1]. The ampoule was embed- 
ded in the same plastic, so that there was no fluid 
between the heating elements and the ampoule. Con- 
vection was measured by taking streak photographs 
of small sulphur particles suspended in organic melts. 
Even though no circumferential variation in tempera- 
ture could be measured, the convection was seldom 
axisymmetric. Furthermore the vertical temperature 
profile in the heater exerted a very strong influence on 
the convection in the melt. If the temperature in- 
creased with height, the convection became immea- 
surably small. If the temperature decreased with 
height, the convection was very vigorous. 

The purpose of the research reported here was to 
determine the influence of heater temperature profile 
on impurity segregation and to relate this to the 
convection in the melt. 

2. EXPERIMENTAL METHODS 

The experimental methods are described in detail 
elsewhere [12]. Naphthalene doped with approx. 
0.3% m azulene was used for the experiments. Azu- 
lene has an equilibrium distribution coefficient of 
about 0.3 in naphthalene [13]. The Schmidt number 
was estimated to be 420 [12]. 

The primary heater was constructed by wrapping 
nichrome wire around a 23 cm long Pyrex tube. An 
additional glass tube was put inside this tube and 
another outside it, all fitting together snugly. Below 
the primary heater was a shorter 6 cm long heater 
constructed in a similar fashion and with a separate 
power supply. Below the heaters were a layer of 
insulation and then a circulating constant tempera- 
ture bath. 

The azulene-doped naphthalene was contained in 
15.6 mm i.d. Pyrex ampoules. The ingot resulting 
from directional solidification was transparent except 
near the top end. This ingot was removed from the 
ampoule by heating until it slid out. It was cut into 
sections with a saw. Samples were taken by scratching 


local areas with glass capillaries, each removing 
about 0.1 mg of material. The samples were dissolved 
into acetone and the amount of naphthalene and 
azulene determined by gas chromatography using a 
flame ionization detector. 

For visualization of convection, a similar appar- 
atus was built with provision for laser slit illumina- 
tion, as in Refs [10] and [11]. Convection was revealed 
by taking time exposure photographs of small sul- 
phur particles suspended in salol (phenyl salicylate). 

3. EXPERIMENTAL RESULTS FOR STABILIZING 
THERMAL CONDITIONS 

By adjusting the temperatures of the two heaters a 
wide variety of temperature profiles were achieved. 
These can be classified into two caregories; those in 
which the temperature increased with height (stabiliz- 
ing) and those in which the temperature decreased 
with height (destabilizing). Figure 1 shows the tem- 
perature profile used for all of the solidification 
experiments performed under stabilizing conditions. 
This profile was measured with thermocouples 
cemented to the outer wall of the Pyrex ampoule. The 
radial variation in melt temperature near the 
solid-liquid interface was about 0.2 K, which corre- 
sponds to a Grashof number of about 10. In flow 
visualization experiments performed under similar 
conditions, significant convection was observed only 
near the top of the ampoule where the temperature 
decreased with height. In the lower portion of the 
ampoule the sulphur particles settled slowly, but their 
paths showed slight deviations from the vertical. 
From these deviations it was estimated that convec- 
tive velocities near the solid-liquid interface were of 
the order of 10^m/s. 

Figure 2 shows the axial azulene concentration 
profile for a freezing rate of about 3 mm/h (0.8 ^m/s), 



the melt position is that at the beginning of a growth run; 
during growth the melt level slowly falls as the ampoule is 
lowered. 
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Fig. 3. Diamteric concentration profiles for the same ingot 
as in Fig. 2, at different axial positions. Here g is the mass 
fraction of the original charge that had been solidified at the 
position. 


which corresponds to a Peclet number Pe of 3.2 
(Pe = RVjDy where R = 7.8 mm and D was estimated 
to be 2 x 10 -5 cm/s). This concentration profile is 
typical of what one would expect to obtain with 
significant convection in the melt during solidifica- 
tion. Figure 3 shows variations in azulene concen- 
tration along a diameter at different positions down 
the ingot of Fig. 2. (Here g is the mass fraction 
solidified.) Because the interface was nearly planar, 
these variations in concentration could only have 
been caused by convection. Since azulene is rejected 
by the growing naphthalene, its concentration in- 


creases as the melt moves along the freezing interface. 
Thus we expect that convection of the melt along the 
freezing interface was in the direction of increasing 
azulene concentration in the resulting solid. In this 
way we can conclude from Fig. 3 that the convection 
was asymmetric during solidification of this material. 

We define two parameters to characterize the cross 
sectional variation in azulene concentration and its 
asymmetry. The “segregation index” is defined as the 
absolute value of the maximum deviation in concen- 
tration from the mean, divided by the mean. Thus in 
Fig. 3 for g = 0.7 the segregation index is about 0.45. 


Table 1 . Results from experiments performed under the thermally stabilizing condition shown 

in Fig. 1 

Pcclet No. 

Crystal 

Mass fraction 

Segregation 

Asymmetry index 

3.2 

30 

0.2 

0.27 

0.51 



0.4 

0.19 

0.34 



0.7 

0.45 

1.56 


32* 

0.2 

0.44 

0.84 



0.6 

0.11 

0.17 



0.7 

0.68 

2.8 


53 

0.2 

0.13 

0.13 



0.5 

0.31 

0.50 




0.17 

0.02 

5.6 

35 

0.6 

0.22 

0.05 


54 

0.5 

0.30 

0.74 




0.49 

0.76 


37t 

0.3 

0.64 

1! 


0.6 

0.42 

0.86 

7.4 

50 

0.2 

0.46 

0.44 




0.29 

0.20 



0.5 

0.57 

0.89 




0.45 

0.17 


55 

0.5 

0.59 

0.71 



0.59 

0.45 

8.6 

34 

0.4 

0.45 

0.21 


0.6 

0.21 

0.02 


49 

0.2 

0.29 

0.06 



0.28 

0.20 



0.5 

0.67 

1.3 




0.51 

0.13 


•Ampoule rotated at 2/3 rpm during solidification, 
tlnsulation applied to one side of apparatus to yield a temperature difference or 6.4 K across 
the outside of the ampoule near the solid-liquid interface. 
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The “asymmetry index” is defined as the maximum 
ratio of concentration at + R to concentration at 
— R, minus 1. In Fig. 3 the asymmetry index for 
g — 0.7 is about 1.56. The segregation index would be 
0 for a constant concentration, while the asymmetry 
index would be 0 for a completely symmetrical con- 
centration profile. 

Table 1 summarizes the results for thermal stabiliz- 
ing conditions. No trends are obvious above the large 
scatter in values of segregation and asymmetry indi- 
ces. This reflects the stochastic nature of the convec- 
tion pattern produced in each experiment. Note that 
segregation and asymmetry changed dramatically 
down each ingot, reflecting the change in convection 
pattern during solidification. 

One would expect that at sufficiently high freezing 
rates the mass transfer in the melt would become 
diffusion controlled, causing the segregation and 
asymmetry indices to decline and approach the values 
expected for diffusion controlled conditions. 
Although the data suggest there may be a maximum 
in segregation index at a Peclet number of about 7, 
there is too much scatter to be certain. 

Thermocouple measurements around the outer 
wall of the ampoule in these experiments showed no 
circumferential variation in temperature, within a 
sensitivity of 0.1 K. It is commonly believed that a 
constant temperature circumferentially must yield 
axisymmetric convection and impurity concentration 
patterns. Therefore an experiment was performed in 
which the ampoule was rotated at 2/3 rpm, in order 
to further reduce any circumferential variation in 
temperature. As shown in Table 1, the segregation 
and asymmetry indices were unchanged, within ex- 
perimental scatter (which was large). Similarly, insu- 
lating one side of the heater to deliberately produce 
an asymmetric temperature field did not significantly 
increase the segregation and asymmetry indices. 



FRACTION SOLIDIFIED 

Fig. 4. Axial concentration profiles for ingots solidified 
at three different rates under the thermally stabilizing 
condition. 



Fig. 5. Typical axial temperature profile for thermally 
destabilizing condition. Based on the temperature difference 
between the maximum and the top of the melt, the Grashof 
number is estimated to have been 158,000. 


Axial concentration profiles also reflect the inten- 
sity of convection. Figure 4 shows such profiles for 
three different freezing rates. The concentration 
profile at Pe = 3.2 is characteristic of convection- 
controlled mass transfer. As Pe was increased, the 
azulene concentration in the first 80% of the ingots 
increased significantly, indicating the decreasing 
influence of convection. 

4. EXPERIMENTAL RESULTS FOR DESTABILIZING 
THERMAL CONDITIONS 

Figure 5 shows a typical destabilizing temperature 
profile. This profile was measured by a fine thermo- 
couple inserted into the melt. Convection was vigor- 
ous throughout the ampoule, and was rarely 



Fig. 6. Concentration profiles along two perpendicular di- 
ameters at the position where 20% of the melt had been 
solidified. Ingot 51 with destabilizing thermal conditions, 
Gr = 76,200 and Pe = 7.4. 
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Table 2. Results from experiments performed under thermally destabilizing conditions 


Peclet No. 

Crystal 

Grashof No. 

Mass fraction 

Segregation 

Asymmetry index 

3.2 

41 

55,600 

0.5 

0.05 

0.05 


42 

88,500 

0.2 

0.93 

0.07 





0.93 

0.22 




0.5 

0.04 

0.02 





0.04 

0.03 


43 

116,000 

0.5 

004 

0.04 

5.6 

44 

39,600 

0.5 

0.17 

0.36 


45 

121,000 

0.5 

0.07 

0.08 





0.07 

0.07 

7.4 

51 

76,200 

0.2 

0.45 

0.13 





0.25 

0.27 




0.5 

0.14 

0.18 





0.26 

0.68 


52 

138,000 

0.5 

0.13 

0.13 





0.12 

0.07 

8.6 

47 

76,700 

0.5 

0.21 

0.46 


48 

109,000 

0.2 

0.47 

1.2 




0.5 

0.50 

0.61 





0.50 

0.38 


46 

158,000 

0.5 

0.15 

0.35 


axisymmetric. The convection velocities were on the 
order of ! mm/s, or about 100 times larger than with 
a stabilizing profile. The convection increased as the 
temperature gradient was increased. The temperature 
gradient was characterized by a Grashof number that 
used the temperature difference between the tempera- 
ture in the melt and the temperature at the upper 
surface of the melt measured when the top of the melt 
was 10.5 cm above the insulated zone, corresponding 
to g = 0.2. 

Figures 6 and 7 show typical cross sectional vari- 
ations in azulene concentration. Visually the dia- 
metric concentration profiles often appeared nearly 
symmetrical, sometimes with a maximum near the 
center and other times with a minimum. As in the 
stabilizing configuration, the concentration profiles 
changed dramatically down each ingot. 

As shown in Table 2, the segregation and asymme- 
try indices again showed a great deal of scatter, but 
tended to be smaller than with thermally stabilizing 



Fig. 7. Concentration profiles along two perpendicular di- 
ameters in the same ingot as Fig. 6, except at the position 
where 50% of the melt had been solidified. 


conditions. The asymmetry appeared to be smallest 
when Pe was small and Gr was large (low growth rate 
with vigorous convection). 


5. IMPLICATIONS FOR SOLIDIFICATION IN SPACE 

A common justification for use of microgravity in 
directional solidification is that more uniform doping 
levels should be obtained. However the present re- 
sults show that under some conditions inhomogenei- 
ties could be greater in space than on Earth. In 
microgravity, weak buoyancy-generated convection 
will occur. If this weak convection is combined with 
low freezing rates, greatly increased inhomogeneities 
could result. Thus when low freezing rates are necess- 
ary in order to avoid constitutional supercooling, it 
may be necessary to combine microgravity and mag- 
netic fields in order to achieve homogeneous crystals. 
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